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ABSTRACT 

In recent years evidence has accumulated suggesting that the gas in galaxy clusters 
is heated by non-gravitational processes. Here we calculate the heating rates required 
to maintain a physically motived mass flow rate, in a sample of seven galaxy clusters. 
We employ the spectroscopic mass deposition rates as an observational input along 
with temperature and density data for each cluster. On energetic grounds we find that 
thermal conduction could provide the necessary heating for A2199, Perseus, A1795 
and A478. However, the suppression factor, of the clasical Spitzer value, is a different 
function of radius for each cluster. Based on the observations of plasma bubbles we also 
calculate the duty cycles for each AGN, in the absence of thermal conduction, which 
can provide the required energy input. With the exception of Hydra- A it appears 
that each of the other AGNs in our sample require duty cycles of roughly 10 6 — 
10 7 yrs to provide their steady-state heating requirements. If these duty cycles are 
unrealistic, this may imply that many galaxy clusters must be heated by very powerful 
Hydra-A type events interspersed between more frequent smaller-scale outbursts. The 
suppression factors for the thermal conductivity required for combined heating by 
AGN and thermal conduction are generally acceptable. However, these suppression 
factors still require 'fine-tuning' of the thermal conductivity as a function of radius. 
As a consequence of this work we present the AGN duty cycle as a cooling flow 
diagnostic. 

Key words: hydrodynamics, cooling flows, galaxies: active, galax- 
ies:clusters:individual: (Virgo, A2199, Perseus, Hydra-A, A2597, A1795, A478) 



1 INTRODUCTION 

The radiative cooling times of the hot, X-ray emitting gas 
in the centres of many galaxy clusters may be as sho rt as 
10 6 years. In the classical model (see e.g. Fabian 1994, for a 
review), where radiative losses occur unopposed, a cooling 
flow develops in which the gas cools below X-ray temper- 
atures and accretes onto the central cluster galaxy where 
it accumulates in molecular clouds and subsequently forms 
stars. However, studies have demonstrated both a lack of 
the expected co ld gas at temperatures well below l-2keV 
(e.g. Edge 2001) and that spectoscopically determined mass 
deposition rates are a f actor of 5-10 less tha n the classi- 
cally determined values JVoigt fc Fabianll2004lh These find- 
ings suggest that the cluster gas is being reheated in order 
to produce the observed minimum temperatures and the re- 
duction in excess star formation compared to expectations. 
As yet it is unclear whether this heating occurs continuously 
or periodically. 

* E-mail:edpope@soton. ac.uk 



There is observational evidence that the temperature 
profiles of galaxy cluster atmospheres are well described by 
the same mathematica l function across a range in redshift 
(e.g. lAllen et alJ 120011) . Given such similarity it is possible 
that the atmospheres of galaxy clusters are in a quasi-steady 
state and that observable parameters such as the radial tem- 
perature profile do not vary signicantly over the lifetime of 
a cluster. This universal temperature profile would be diffi- 
cult to sustain if the density profiles varied significantly with 
time. In addition, if large quantities of gas were deposited in 
the central regions one might expect that the density pro- 
files would be much more centrally peaked than observations 
suggest. The culmination of this argument is that, at least 
in the central regions, the inward flow rate of mass must be 
roughly independent of radius to ensure that mass is not de- 
posited at any particular location which would significantly 
alter the density profile. 

One particular problem with this is that all of the mass 
flowing in this region must be deposited at the cluster centre. 
Therefore, unless the mass flow rate is small this would result 
in large quantities of gas at the cluster centre. 
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In contrast, spectroscopically determined mass depo- 
sition rates suggests that mass is deposited at a roughly 
constant rate throug hout the cluster, at resolved radii 
llVoigt fe Fabian|l2004h . However, a roughly constant mass 
deposition rate is indicative of a mass flow rate that is 
roughly proportional to radius, rather than constant as the 
density profiles suggest. Yet, if such a flow were persistent 
this proportionality would result in density profiles which 
contradict the observational evidence by being more cen- 
trally peaked. 

If we are to reconcile the implications of the density pro- 
files and mass deposition rates then one possible explanation 
is that the mass flow rate exhibits two different asymptotic 
properties. That is, near the cluster centre or within the cen- 
tral galaxy the mass flow rate should be roughly constant so 
that the density profiles are essentially left unchanged, while 
the mass deposited at the cluster centre is sufficiently small 
to agree with observations. At larger radii the mass flow rate 
should still be proportional to radius and satisfy the obser- 
vations on resolvable scales which imply a constant mass de- 
position rate. The relationship between the mass flow rates 
and the mass deposition rates are defined section 3. 

The two main candidates for heating clus- 
ter at mospheres are: Active Galactic Nuclei (AGNJ 
fe.g. iTabor fc Binnevl Il993t IChurazov et all l200ll 
iBriiggen fc Kaiser! I200S : iBriiggenl 12003ft and thermal 
conduction (e.g'H GagtJ f98^T| za]^^ka fc: Naravanl l2003t 
IVoigt et al-lboollVoigt fc Fabianl \2004) . Heating bv AGN 
is thought to occur through the dissipation of the internal 
energy of plasma bubbles inflated by the AGN at the centre 
of the cooling flow. Since these bubbles are less dense than 
the ambient gas, they are buoyant and rise through the 
intracluster medium (ICM) stirr ing and exciting sound 
waves in the surounding gas (e.g. Ruszkow ski fc Begelmar] 
120021 : iFabian et al.ll2002(i . This energy may be dissipated by 
means of a turbulent cascade, or viscous processes if they 
are significant. Deep in the central galaxy other processes 
such as supernovae and stellar winds will also have some 
impact on the ambient gas. 

However, AGN are only periodically active which could 
result in similarly periodic heating rates, thus making the 
possibility of a totally steady state unlikely. In this case 
one could imagine a scenario in which a quasi steady-state 
is possible where temperature and density profiles oscillate 
around their average values. It is also possible that dissipa- 
tion of the plasma bubbles may occur over timescales longer 
than the A GN duty cycle thus pr oviding almost continous 
heating (e.g. lRevnolds et al]l2005h . 

Thermal conduction may also play a significant role in 
transferring energy towards central regions of galaxy clusters 
given the large temperature gradients whi ch are observed in 
many clusters. In fact , several authors (e.g. lVoigt et all2002t 
IVoigt fc FabiaT] 12004 have shown that on energy grounds 
alone it may be possible for thermal conduction to provide 
the necessary heating in some clusters. 

However, the exact value of the thermal conductivity re- 
mains uncertain. The theoretical value for the thermal con- 
ductivit y of a fu l ly ion ised, unmagnetised plasma is calcu- 
lated bv ISpitzerl (I1962T) . Magnetic fields, which are thought 
to exist in the cluste r gas, based on Farad ay rotation mea- 
sures of clusters (e.g. lCarilli fc TavloJl20 02'l . may greatly al- 
ter this value. The standard method by which the unknown 



effects of magnetic fields are taken into account is to include 
a suppression factor, /, in the Spitzer formula, which indi- 
cates the actual value as a fraction of the full Spitzer value. 

Several authors have found temperature profiles in some 
galaxy clusters that are compatible with thermal conduc- 
tion, and values of t he suppression factor that a re physi- 
cally meaningful (e.g. IZakamska fc Naravanl I2003F) . In con- 
trast, more detailed work in which the Virgo cluster was 
simulated using complete hydrodynamics, including radia- 
tive cooling and heating by thermal conduction has demon- 
strated that thermal conduction cannot prevent the oc- 
cur ence of a cooling catastrophe, in this particular exam- 
ple iPope et alJl20031 . Furthermore, the mass flow rates un- 
der such circumstances are not constant with radius so that 
mass builds up in particular areas and therefore changes the 
density distribution, making it much more centrally peaked 
than observed in real clusters. 

To answer the questions regarding the necessary heat- 
ing rates and mechanisms involved in galaxy clusters we 
construct a simple model based on two main assumptions: 
galaxy clusters are in approximate steady state, and the 
mass flow rates fulfill the constraints provided by observa- 
tions at both large and small radii. To ensure this, the radial 
behaviour of the mass flow rates is modelled using an ob- 
servationally motivated, but empirical function consistent 
with the observations. Using mathematical functions fitted 
to the temperature and density data of seven galaxy clusters 
we derive the required heating rates within regions which are 
consistent with those for which the mass flow rates were de- 
termined. We also compare the observed energy currently 
available in the form of plasma bubbles with the heating 
requirements for each cluster. From the heating rate pro- 
files we calculate the thermal conduction suppression factor 
as a radial function for each cluster in order that we may 
determine which clusters could be heated by this process. 

The plan for this paper is as follows. In Section 2 we 
give the details of the parameters used to fit observational 
temperature and density data for seven galaxy clusters. In 
Section 3 and 4 we describe the model used to estimate the 
required heating, and thermal conduction suppression fac- 
tors. The heating rates we have calculated are compared 
with observational estimates of both AGN heating and a 
combination of AGN heating and thermal conduction in Sec- 
tions 5 and 6. In Section 7 we summarize our main findings. 

The results are given for a cosmology with Ho — 
TOkms^Mpc" 1 , f2 M = 0.3 and fi A = 0.7. 



2 FUNCTIONS FITTED TO OBSERVATIONAL 
DATA 

In this section we give details of the functions used to derive 
the mass deposition and heating rate profiles, the references 
to the observational data, and show a comparison with the 
data in figure Q 

We have only fitted our own functions to the data when 
the original authors had not included confidence intervals for 
their model parameters, but we give the functions used to 
fit all of the cluster atmospheres. The values of the fitted 
parameters are given in tables and [5] 
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Name 


no (cm 3 ) ni(cm 3 ) 


A) 


Pi 


ro(kpc) 


n (kpc) 


Virgo 


0.089±0.011 0.019 ±0.002 


1.52±0.32 


0.705 


5±1 


23.3±4.3 


Perseus 


0.071±0.003 


0.81±0.04 




28.5±2.7 




Hydra 


0.07±0.02 


0.72±0.004 




18.6±0.5 




A2597 


0.13±0.07 


0.15±0.10 


0.79±0.06 


1.0±2.6 


43±14 


A2199 


0.19 


0.75 




1 




A1795 


0.066±0.067 


0.24±0.11 


0.41±0.13 


57±37 


12±4.9 


A478 


0.153±0.019 


0.55±0.05 


0.41±0.34 


145±32 


6.6±1.9 



Table 1. Summary of best fit parameters for density profiles with 1-sigma errors obtained from the least- 
squares fitting procedure. We do not give errors the model function parameters for A2199 since we were 
unable to obtain the observational data. See table 2 for references. 



The Virgo cluster's electro n number density data was 
fitted by Ghi zzardi et alJ J2004I1 using a double ,3-profile, 



no 



- + (1) 

Fo r simplicity, the Pe rseus (see lSanders et al.l2004l) and 
Hydra JPavid et alJl200lfl density data were fitted with sin- 
gle /3-profiles, 



no 



[1 + (r/r ) 2 ]^ • 



(2) 



We were unable to obtain the observational data for 
A2199 a nd so relied on the fit s presented by the original 
authors jjohnstone et all I2002T) who found that a simple 
power-law was sufficient to describe the density distribution, 



no 



-00 



(3) 



For A2597, A1795 and A478 (s ee lMcNamara et all200ll 
lEttori et alj|2002t ISun et alJl200l respect i vely^ i we employ 
the function used bv TPemrisfcCThandranl i2005l) to fit the 
density distributions, 



no 



[1 + {r/roYYo [1 + (r/n) 2 ]/3i 



(4) 



As above, for the cluster gas temperature profiles we 
have only fitted functions to the data when the original au- 
thors had not included confidence intervals, but give the 
det ails for all cluster s in ou r sample. 

iGhizzardi et al.l <l2004h found that the Virgo tempera- 
ture data is best described by a Gaussian curve, 



T = T 



1 r 2 



(5) 



IChurazov et alJ <l200ot) fitted the Perseus temperature 
data with the function, 



We fitted the temperature data for A2597, A1795 
and A 478 using the same function as iDennis fc Chandranl 
(2005). 



T = To 



[1 + (r/r ct ) 2 } 



21 S ■ 



(T) 



The b est-fit tempe r ature profile for the Hydra cluster 
is given bv lDavid et alJ l|200lF ) as a power-law, 



T = T 



(8) 



The A2199 temperature d ata was found by the original 
authors jjohnstone et al]l2002T) to be described by a power- 
law without a characteristic length-scale. For consistency we 
use equation JSJ to represent the temperature distribution 
in this cluster. 

Although we have provided 1-sigma errors for the fit- 
ted parameters that describe the temperature and density 
distribution of the gas in galaxy clusters, we will not give 
errors for any quantities derived from these parameters. The 
reason for this is that the errors and gradients of the tem- 
perature and density are all model dependent. In addition, 
the calculations are complex, making the standard propa- 
gation of errors impractical and error estimates themselves 
very uncertain. In addition, although we can make very con- 
servative estimates of the confidence limits, of a particular 
derived quantity, they are too vague to be of importance. 



3 THE MODEL 

3.1 General heating rates 

Starting from the assumption that the atmospheres of 
galaxy clusters are spherically symmetric, and in a quasi 
steady-state, it is possible to derive what the time averaged 
heating rate, as a function of radius, must be in order to 
maintain the observed temperature and density profiles, 



T = T 



l + (r/r ct y 



5+(r/r ct y 



(6) 



2 A 1 d 

h = n A rad + ^73- 



M ( 5k b T 
47T y 2/im p 



(9) 
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Name T (keV) Ti(keV) r c t(kpc) <5 Reference 



Virgo 


2.4±0.1 


0.77±0.10 


23±4 




1 


Perseus 


8.7±1.1 




63±5 


2.7±0.3 


2 


Hydra 


2.73±0.07 




10 


0.12±0.01 


3 


A2597 


4.1±0.1 


2.0±0.2 


16±10 


0.7±0.4 


4,5 


A2199 


1.0 




1 


0.29 


6 


A1795 


7.2±1.9 


5.1±2.2 


40±34 


0.5±0.7 


7,5 


A478 


9.7±1.5 


7.6±1.7 


16±5 


0.3±0.2 


8,5 



Table 2. Summary of temperature parameters fitted to the data with 1-sigma errors obtained from the 
least-squares fitting procedure. References:- (1) Ghizzardi et al. (2004); (2) Sanders et al. (2004); (3) David 
ct al. (2001); (4) McNamara et al. (2001) and RafFerty et al. (2005), in preparation; (5) Dennis & Chandran 
(2005); (6) Johnstone et al. (2002); (7) Ettori ct al. (2002); (8) Sun et al. (2003). We do not give errors for 
A2199 since we were unable to obtain the observational data. We also do not give an error for the Hydra 
temperature scale- height since none was given in the original publication, David et al. (2001). 




1 10 100 1 10 100 

distance [kpc] distance [kpc] 



Figure 1. Left: density data for all clusters except A2199 for which we show the best-fit, since we did not have access to the observational 
data. Right: temperature data except for A2199 and Hydra for which we give the fitted functions, since we did not have access to the 
observational data. The key in the right-hand panel applies to both temperature and density data. 



where M is the mass flow rate, T is the gas temperature, n 
is the gas number density, A ra( j is the cooling function and <f> 
is the gravitational potential. Here we assume only subsonic 
flow of gas. 

The mass flow rate through a spherical surface, radius 

r, is 

M = Airr 2 pv T , (10) 

where p is the gas density and v T is the gas velocity in the 
radial direction. 

The difference between the rate at which mass enters 
and leaves a spherical shell, of thickness Ar, is called the 
mass deposition rate, AM. In terms of the mass flow rate, 
the mass deposition rate is 

AM«(f)*r (ID 



Consequently, the mass flow rate, at a particular radius, 
can be recovered by summing the contributions of the mass 
deposition rate in each shell up to that radius. Henceforth, 
we refer to the mass flow rate as the integrated mass depo- 
sition rate. 

For subsonic gas flow the cluster atmosphere can be 
assumed to be in approximate hydrostatic equilibrium. On 
this basis, we are able to estimate the gravitational accelera- 
tion from the temperature and density profiles of the cluster 
atmosphere, 



dP _ _ d0 
dr ~ P dr' 

where P is the gas pressure and p is the gas density. 
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The gas pressure is related to density and temperature 
by the ideal gas equation 



Therefore, the rate at which mass is deposited within 
the jth shell is 



P = 



pk h T 
Atm p ' 



(13) 



where (jm p is the mean mass per particle. We assume that 
/i=16/27 as appropriate for standard primordial abundance. 

To avoid anomalies when calculating spatial derivatives, 
we fit continuous analytical functions through the density 
and temperature data in the previous section. This ensures 
that we do not encounter any large discontinuities which 
may result in extreme heating rates. 

The only unknowns in equation © are then the inte- 
grated mass deposition rate, M, and the heating rate, h. 



3.2 Calculating Mass Deposition rates 

In the classical cooling flow picture, which assumes that the 
heating is zero everywhere, and ignoring the effect of the 
gravitational potential, the bolometric X-ray luminosity, in- 
tegrated mass deposition rate and gas temperature are re- 
lated by integrating equation (e.g. lFabianlll994f) over a 
spherical volume. Analytically, the bolometric X-ray lumi- 
nosity within a given radius is simply the volume integral of 
the radiative cooling rate per unit volume, 



Lx(< r) = 47r 



n 2 A rad r 2 dr. 



(14) 



Integrating the second term in equation @ gives the 
reationship between bolometric luminosity, integrated clas- 
sical mass deposition rate and temperature, 



Lx(<r) = 



5fc b 
2fim p 



M c i aa (r)T(r). 



(15) 



For observational data the luminosity and classical mass 
deposition rates at each radius are estimated by dividing 
the observed 2-d projection of the cluster into concentric 
shells. Using this method one is able to reconstruct the 3-d 
properties of a cluster. The contribution, to the total X-ray 
luminosity, of the jth spherical shell, of thickness Ar, is then 
calculated by discretizing equation 1151 . 



5fcb 
2/im p 



(MdasjATj + AM clasj Tj) 



(16) 



In any given shell, the term involving AT) is assumed 
to be small compared to the term involving AAf c ; as j . By ig- 
noring the ATj term the resulting mass deposition rate is 
larger than it would be if the change in temperature across 
the shell was taken into account. Therefore, the classical 
value is the absolute maximum possible mass deposition 
rate. T he luminosity of a gi ven shell is then assumed to 
be fe.g.lVoirt Fahi^EifMh. 



AL 



Xj 



; 4nr 2 Arn 2 A ra( i 



5fc b Tj 
2(im p 



AM, 



clasj ■ 



(17) 



AM, 



clasj — 



47rr-j 2 Arnj 2 A r 

2/j,m,p 



(18) 



The total X-ray luminosity emitted by the gas within a 
particular radius is calculated by summing the contributions 
from all of the shells within this radius 



Lx(< r) 



5fcb 
2/im p 



M c i as (<r)r(r), 



(19) 



where Af c i as (< r) is the integrated mass deposition rate ob- 
tained by summing the contribution from each shell within 
the radius, r. 

If we include a heating term the expression equivalent 
to equation (1151 becomes, 



Lx(< r)-H(<r) 



5fc b 
2(im p 



M coo iT(r), 



(20) 



where M coo i is the integrated mass deposition rate, resulting 
from the excess of cooling over heating. This integrated mass 
deposition rate is distinct from the integrated classical mass 
deposition rate which is only used to describe the situation 
in the absence of heating. 

With equation 1161 becomes, 



A L 



5fcb 
2/im p 



(Mcooij ATj + AAfcooijTj) . (21) 



In practice, a realistic estimate of the mass deposited 
into each shell, which takes into account heating, can be 
obtained by fitting models to the X-ray spectrum of each 
shell or a larger region if desired. The spectroscopic mass 
deposition rate is just a fitting parameter in this case and 
an integrated value can be obtained in exactly the same way 
as for the classical case above. 

Generalising equation ll9H to allow for heating we find 
that if the spectroscopically determined mass deposition 
rate, M b s , is representative of the actual mass deposition 
rate, M coo i, then one can determine the heating rate from, 



L x (< r)-H(<r) = 



5fcb 
2/im p 



M obs (< r)T(r) = — — M cooI T(r) 

(22) 

If the spectroscopically determined value for the mass 
deposition rate is a true representative of the real mass de- 
position rate, then the integral within the cooling radius is 
simply the integrated mass deposition rate within the cool- 
ing radius, irrespective of the radial distribution. Therefore, 
it would seem a sensible starting point for any calculations 
of the heating rates in galaxy clusters. To ensure consistency 
between theory and observations, we must employ the same 
mass flow rate (integrated mass deposition rate) over the 
same region, in our heating models. Consequently, if we use 
the observed spectroscopic mass deposition rate integrated 
up to the cooling radius, then to be strictly accurate, the 
heating rates we subsequently calculate are only valid up to 
this radius. 
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3.3 Model Integrated Mass Deposition Rates 



3.4 Thermal conduction 



We now construct a model for the integrated mass deposi- 
tion rate in galaxy clusters on the basis of the observational 
evidence. Although the model function is essentially empiri- 
cal there are three main constraints to which it must adhere. 
These are: firstly, the fact that the observed spectroscopic 
mass deposition rate in each shell is approximately constant 
up to the cooling radius suggests linear dependence of flow 
rate on radius. Secondly, the observed density profiles also 
suggest that the mass deposition rate is negligible near the 
cluster centre. Otherwise a pronounced peak would be ob- 
served in the density profiles. Thirdly, the integrated mass 
deposition rate at the cooling radius must be equal to the 
observed integrated mass deposition rate at that location. 
A simple empirical form for the integrated mass deposition 
rate which is, in some way, consistent with the above is 



M(r) = M obs {< r cool ) 



K + 



rcool 



Thermal conduction of heat from the cluster outskirts 
to their centres may provide the required heating of the 
central regions without an additional energy source, like 
an AGN (e.g. iGaetd Il989l: IZakamska fc Naravanl 120031: 
IVoiet fc Fabianll2004h . Of course, thermal conduction only 
transports energy from one region to another and does not 
lead to net heating. However, if we consider the case for 
which there is an infinite heat bath at large radii, then any 
drop in temperature at large radii, due to the inward trans- 
fer of thermal energy, is negligible. 

Several 1-d models have been able to achieve 
a steady-state using ther mal conduction alone (e.g. 
Zaka mska fc Naravanl EOO^ and the combined effects 



(23) 



of thermal conduction and AGNs (e.g. Briiggen 2003; 
iRuszkowski fc Begelmanll2002D . We, therefore, compare the 
required heating rates with an estimate of the energy trans- 
port supplied by thermal conduction given the current state 
of each galaxy cluster. 

The volume heating rate for thermal conduction from 
an infinite heat bath is given by 



1 d 



where M b s (< r coo i) is the integrated spectroscopic mass 
deposition rate at the cooling radius, K is a constant and 
fcooi is the cooling radius. 

To estimate K we express it the following way, K — 
(nt/Vcooi) 2 , where tyl is a length scale corresponding to K. 
A suitable value for tk is the radius of the central galaxy, 
roughly 30 kpc, since the effect of the interstellar medium is 
likely to significantly alter the inflow of material. This leads 
to values for K ranging from 0.73, for Virgo, to 0.04, for 
A478. 

Note that our expression for M implies that AM /dr — 
for samll radii, meaning that no material is deposited. How- 
ever, at or near r — the material must be deposited since 
it cannot flow to negative radii. This can account for the 
excess star formation which is observed to occur in these 
regions. For reasonable gradients of the temperature distri- 
bution, the second term on the right hand side of equation 
© then tends to infinity as r tends to zero. This implies that 
our model predicts negative heating rates, h, for small radii. 
In otherwords, the model is only valid outside the region at 
the very centres of clusters. To indicate the range over which 
our model can provide physically meaningful predictions, we 
list in Table 3, the minimum radius for each cluster at which 
h=0. We take this minimum radius as the lower integration 
limit in our calculation of the total heating rates. 

For each cluster the observations are inevitably limited 
by spatial resolution and in some cases the lower integration 
limits, explained above, occur at radii which are closer to the 
cluster centre than the spatial resolution limits. In principle, 
the spatial resolution presents a limit beyond which we are 
unable to accurately describe the temperature and density 
of the cluster gas. However, because we have fitted analyti- 
cal functions to the temperature and density data for each 
cluster we can extrapolate the behaviour of the tempera- 
ture and density, and hence the heating rates, from the data 
point nearest the cluster centre further towards the centre. 
The results at r < r m i n should be ignored as we have no 
observational information to constrain them. 



£-cond oi I i 

r A dr V dr 



dT 



(24) 



where k is the thermal conductivity and dT/dr is the tem- 
perature gradient. 

W e take the thermal conductivity to be given bv lSnitzeil 
(1962), but include a suppression factor, /(r), defined in 
the introduction, to take into account the effect of magnetic 
fields, 



1.84 x 10" 5 /(r)T 5/2 



In A c 



(25) 



where lnA c is the Coulomb logarithm. 

For pure h ydrogen the Coulomb logarithm is (e.g. 
IChoudhurilll998l) 



A r = 2Ann e 



a 1 \ -3/2 



(26) 



where n c is the electron number density. 

For a steady-state to exist, the heating by thermal con- 
duction must be equal to the heating rate, 



-h, 



(27) 



where h is the required heating rate per unit volume, given 
by equation 

Thi s is essentially the same ene rgy equation as solved 
by (e.g. IZakamska fc Naravanl l2003fl . although we also al- 
low for mass deposition, line-cooling and variations in the 
Coulomb logarithm. We then solve equation H27II . using the 
observed temperature and density profiles, to determine the 
radial dependence of t he suppression factor. This is differ- 
ent to the approach of j Zaka mska fc Naravanl2003i) . In their 
paper they use the same equations to calculate temperature 
and density profiles that are consistent with a constant sup- 
pression factor. These derived profiles are only constrained 
by observations at the minimum and maximum radii acces- 
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sible to observations. We take the density and temperature 
profiles from the observed data and allow the suppression 
factor to vary with radius. 

Our function for the suppression factor is obtained by 
integrating both sides of equation (12711 over a spherical sur- 
face and then rearranging for /, 

f r r 2 h(r)dr 
^ r > = 1.84 x 10- 5 r 2 dT/drT5/2/lnA c ' ^ 

We note that the energy flux due to thermal conduction, 
K&T/Ar, increases with radius so that, for an infinite heat 
bath at large radii, energy must always be deposited at each 
smaller radius, rather than taken away. 



4.3 Required Heating rates and Thermal 
Conduction Suppression factors 

The heating rates are shown in figure |5J for the form of 
integrated mass deposition rates given in equation 1231 . For 
comparison, we also show the volume heating rate profiles 
for unsuppressed thermal conduction. 

Figure shows that the heating rates for the entire 
sample of clusters exhibit similar profiles. This is because 
the radiative losses most strongly depend upon the density, 
therefore so must the required heating rate. Since the ma- 
jority of the density profiles are described by /3-profiles, this 
similarity is expected. 



4 RESULTS: 1 

4.1 Comparison of Luminosities 

To ensure that our estimates are compatible with those cal- 
culated from observations we compare the bolometric X- 
ray luminosities, within the cooling radius, determined from 
our functions fit ted to the data and the observed luminosi- 
ties presented in lBirzan et alJ (|2004h . We assume half solar 
metallic ity for all clusters and use t he cooling function as 
given in ISutherland fc Dopital dl993t) . 

Table|2]demonstrates reasonable agreement between the 
results o btained using the fi tted functions and those pre- 
sented in lBirzan et alJ i2004h . 



4.2 Integrated Mass Deposition rates 

Since the values of the integrated mass deposition rates 
that we use are critical to determining the heating rates, 
we present these first. In table 4 we compare the inte- 
grated spectroscop ic mass deposition rate, M spcc (< r coo i) 
l|Birzan et ail 12004 V with the inte grated cl assical mass de- 
position rate, M c i as (< r coo i) iFabianl|l994l) at the cooling 
radius, and the central integrated mass deposition rate pre- 
dicted by equation H23|l . 

From equation 12011 it is clear that the spectroscopic 
mass deposition rates are indicative of the energy in- 
jected into the cluster. Because of this, the spectroscopic 
mass deposition rate should be less than the classical 
value,depending upon the relative magnitudes of the heat- 
ing and cooling, as is the case. However, because of the time 
taken to dissipate this injected energy into the ambient gas 
this implies that the current spectroscopic mass deposition 
rates are a function of the heating which has taken place 
over the last few times 10 s yrs. 

The central integrated mass deposition rates predicted 
by equation are sufficiently low that the total mass de- 
posited in the cluster centre, over a Gyr, will not significantly 
alter the total central mass. It should be noted that a smaller 
value of the constant, K, in equation l|23fl . would result in 
still smaller central integrated mass deposition rates. 



In comparison, it is clear that the volume heating rates 
for unsuppressed thermal conduction vary from cluster to 
cluster and do not share the same profile as the required 
heating rates. Even for the clusters in which thermal con- 
duction can supply the necessary energy, the different radial 
dependences require that the suppression factor is fine-tuned 
to provide the correct heating rate for all radii. This demon- 
strates the different nature of the physical processes involved 
in radiative cooling and thermal conduction, making a bal- 
ance between the two hard to achieve. 

Using equation 12811 we calculate the thermal conduc- 
tion suppression factor, as a function of radius, for which 
the required heating rate is achieved. These results, plotted 
in figure |3] suggest that thermal conduction could provide 
the required heating rates in A2199, A478, Perseus and very 
nearly A1795. It seems impossible that thermal conduction 
would be able to transport sufficient energy towards the cen- 
tral regions for the remainder of the sample. 

It is evident that the suppression factor profiles are dif- 
ferent in almost every case, although A1795 and A478 ap- 
pear to have certain features in common. The only corn- 
men trait is that the required suppression factor tends to be 
largest near the centre, except for A2597 and Hydra. This 
suggests that even in the regions where the suppression fac- 
tor takes physically meaningful values (/ < 1), in order for 
thermal conduction to provide the necessary heating, the 
parameters which determine the suppression factor require 
unique fine tuning in each cluster. We note that suppression 
factors of greater than unity are unph ysical for lam inar flows 
but are possible in mixing layers (e.g. lCho et"al"ll2003l) . 

Our suppression factors can be compared with those 
found by Zakamska fc NaravarJ a2003Fl who studied four of 
the clusters in our sample. They find suppression factors or 
1.5, 2.4, 0.4 and 0.2 for Hydra, A2597, A2199 and A1795 
respectively. These values agree with our results in that we 
also find that thermal conduction must take unphysically 
large values for Hydra and A2597 and realistic value for the 
remaining two clusters. IVoigt fc Fabianl i2004h also provide 
a suppression factor, roughly 0.2-0.3, for A478 which is con- 
sistent with, albeit lower than,the possible values we have 
found for the same cluster. 
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1 10 100 

distance [kpc] 

Figure 2. Comparison of cooling and heating rate profiles for each cluster. The heating rates fall to zero at a radius which depends 
upon the mass deposition rate. At the point at which this occurs, we also terminate each of the other curves since wc cannot determine 
the gas properties within this radius. 
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Name Lx (^cool 

)/10 42 ergs- 1 Lx(r C ool)/10 42 ergs '"cool 

/kpc r min /kpc 



Virgo 


10.7 


9.8l° ;? 


35 


1 


Perseus 


550 


670+« 


102 


3 


Hydra 


243 


350±" 


100 


3 


A2597 


470 


430l« 


129 


4 


A2199 


156 




113 


1 


A1795 


493 


490±|° 


137 


1 


A478 


1494 


1220l 6 6 ° 


150 


3 



Table 3. Comparison of bolometric luminosities derived from our functions fitted to the data (column 2) 
with those given in Birzan et al. (2004) (column 3), the cooling radius (column 4) and the lower integration 
limits (column 5). 



Name 


A4 P oc(< r cool ) 


*^clas(< r cool) 


-^spcc(< r cool) r J^ ol 


Virgo 


1 S+ 1 - 2 


10 


1.3 


Perseus 


54 +48. 


183 


4.7 


Hydra 


Ml? 


315 


1.3 


A2597 


59 +40 
oa -40 


480 


3.2 


A2199 


2+ 7 
z -1.9 


150 


0.1 


A1795 




478 


0.9 


A478 


150i 6 6 ° 


570 


6 



Table 4. Mass deposition rate parameters (all quantities in solar masses per year). The integrated spectral 
mass deposition rates are taken from Birzan et al. (2004) and the classically determined values are from 
Fabian (1994) and references therein. The central mass flow rates calculated using equation 1231 are also 
shown for comparison. 



5 COMPARISON OF REQUIRED HEATING 
RATES WITH OBSERVATIONS 

5.1 Comparison with AGN 

The time-averaged mechanical luminosity of a central AGN 
in each cluster is estimated using 

aPV 

(L mcch ) = (29) 

where cv=16. P is the ambient cluster gas pressure, V is 
the bubble volume and t is the estimated duty cycle for 
an AGN. We assume initially that the duty cycle of each 
AGN is 10 s yrs. The factor 16 arises from 2 bubbles per 
outburst, a factor of 2 for the energy dissipated in the shock 
expansion of the bubble and 4 for 7/(7 — 1) where 7 = 4/3 
for relativistic gas. 

This is simply an estimate of the rate at which energy 
is injected into the cluster, by the central AGN, and is not 
related to any particular physical process by which this en- 
ergy is dissipated e.g. the viscous dissipation of sound waves 
or effervescent heating. 



The estimated AGN power injection rates are given in 
column 2 of Table 5. To compare the required heating rates 
with the observational estimates of AGN heating in equa- 
tion 1291 1 we must calculate the volume integral of equation 
@. From this we can estimate the heating luminosity re- 
quired within a given volume to satisfy the requirements of 
the assumed steady-state. Comparing this value with that of 
the power output from AGNs at the centres of galaxy clus- 
ters allows us to estimate the frequency with which buoyant 
plasma bubbles, identical to those currently observed, must 
be produced in order to provide sufficient heating within the 
cooling radius. This is defined as the AGN duty cycle, which 
may be calculated using 

4tt p 001 r 2 hdr t 

7r S = ~. (30) 

where r coo i, the upper integration limit, is the cooling radius 
and r m i n is the radius at which the heating rate tends to zero 
for our chosen descriptions of the mass flow rate, r is the 
calculated duty cycle and t is the time defined in equation 
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Figure 3. Thermal conduction suppression factors for the sample of seven clusters. The thick, solid line shows the maximum physically 
meaningful suppression factor (f=l) for comparison. 



5.2 The Combined Heating by Thermal 
Conduction and AGN 

To estimate the total rate of energy injection by AGNs into 
the cluster gas within the cooling radius we assume that all 
of the energy available in the bubbles is dissipated within the 
cooling radius. Furthermore, the integral of equation 1241 
over a spherical surface gives the rate at which thermal con- 
duction transfers heat across a radius, r, 



L cond (r) = 4nr 2 k,— . (31) 
dr 

The maximum total rate at which energy can be in- 
jected into this region is the sum of these two processes. We 
compare this sum with the heating rate required for steady- 
state and allowing for mass deposition. Within the cooling 
radius, we can then determine the thermal conduction sup- 
pression factor required to maintain a steady-state given the 
current AGN heating rate. This suppression factor is calcu- 
lated using 



(< r ) - (£mech) 



where Lhcat(< T ) is the required heating rate for the spheri- 
cal volume within r, (L mec h) is the current mechanical power 
of the AGN assuming a duty cycle of 10 8 yrs. 



6 RESULTS: 2 

6.1 AGN Duty Cycles 

The duty cycles (see Table 5) for Virgo and A478 are of 
the order of 10 6 yrs which is very short compared to the 
predicted lifetimes of AGN (e.g. iNipoti fc Binnevl2005l and 
references therein). In contrast, the Hydra cluster requires 
recurrent outbursts of magnitude similar to the currently 
observed one only every 10 s yrs, or so. The required duty 
cycles for the remaining AGNs are of the order of 10 7 yrs. 
From this we conclude that if thermal conduction is negligi- 
ble and if this sample is representative of galaxy clusters in 
general, then many if not all clusters will probably be heated, 
at certain points in time, by extremely powerful AGN out- 
bursts. However, it is worthwhile pointing out that roughly 
71% of cD g alaxies at the centres of clusters are radio-loud 
<BurnJl99(l) which is larger than for galaxies not at the cen- 
tres of clusters. This may suggest that the galaxies at the 
centres of clusters are indeed active more frequently than 
other galaxies. It is possible that a combination of these two 
effects satisfies the heating requirements we have identified. 

In contrast to the spectroscopically derived mass de- 
position rates, the plasma bubbles, that still exhibit radio 
emission and are currently observable, are probably more 
indicative of the current heating rate rather than the heat- 
ing rate over the previous 10 9 years, or so. If this is true, 
only Hydra is currently being excessively heated in our small 
sample. 

Our results also show that the gas outside the cool- 
ing radius requires significant amounts of energy to allow a 
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steady-state. Given the large heating requirements for the 
regions outside the cooling radii of Virgo, A2199 and A478, 
it is reasonable to ask whether it is possible for the observed 
plasma bubbles in these clusters to deliver e nergy at such 
large radii. The observational evidence (e.g. IfBirzan et alJ 
|2004) seems to show that bubbles rarely rise to beyond 40 
kpc from the cluster centre and in general tend to reach 
only a few tens of kiloparsecs. This may be because be- 
yond these radii, hydrodynamic inst a bilities have shredded 
the bubbles (e.g. |Kaiser_et_al ]|2005t iRevnolds et al]l2005l 
iRobinson et al.l2004l : IChurazov et al]l200ll) . or alternatively 
the isoentropy radius at which the bubbles spread out and 
form pancakes is gene rally at radii of these magnitudes (e.g. 
IChurazov et al]l200lf) . Alternatively, the apparent absence 
of bubbles at larger radii may be a selection effect since they 
are harder to detect against the rapidly declining X-ray sur- 
face brightness profiles. In any case, it is possible that ex- 
tremely powerful outbursts with highly extended jets such 
as Hydra-A or Cygnus-A type events are the best mecha- 
nism by which energy is supplied to such large radii, since 
they cannot be reached by smaller scale outbursts. 

6.2 The Combined Heating by Thermal 
Conduction and AGNs 

From table|5]it is clear that A2597 is currently being heated 
insufficiently by the combined effects of AGN and thermal 
conduction. Thus, it seems that the only way that A2597 can 
avoid a cooling catastrophe is by an exceptionally large AGN 
outburst which will have to be roughly 100 times more pow- 
erful than the current outburst in order to provide sufficient 
heating. Of the other clusters in this sample, the suppres- 
sion factors are all physically acceptable except for Hydra. 
The reason for this and the negative value of the suppression 
factor for Hydra is that the energy of the current AGN out- 
burst exceeds the heating requirements within the cooling 
radius. However, the errors associated are relatively large 
so it is difficult to say with any certainty if a combination 
of heating by thermal conduction and AGN can achieve a 
steady-state. 

Although the global suppression factor for thermal con- 
duction is physically acceptable, it still requires fine-tuning 
inside each cluster to allow a steady-state (see Section 4.3). 



7 SUMMARY 

In the heating of intracluster gas, the energy input required 
from AGNs is likely to depend on the strength of any other 
heating processes. For example, if thermal conduction is 
present and can reduce the integrated mass deposition rate 
towards the cluster centre, it is likely that in this case the 
AGN outbursts will have to be less powerful than in an oth- 
erwise identical cluster in which thermal conduction does 
not occur. There may also be additional heating processes 
such as stirring by galaxy motions or the gas motion result- 
ing from past cluster-cluster mergers. 

We have calculated the heating rates required to main- 
tain a steady state for a sample of seven galaxy clusters with 
the assumption that the mass flow rates are independent of 
radius. For our model we use the spectroscopically deter- 
mined integrated mass deposition rates. Here we summarise 



our main findings for each cluster in our sample in terms of 
the mass deposition rates, thermal conduction and required 
AGN duty cycles. 

A2199: It appears that, without AGN heating, thermal 
conduction is sufficient for the steady-state heating require- 
ments. In the absence of thermal conduction, the required 
AGN duty cycle of roughly 10 7 yrs is short, but still com- 
patible with predicted AGN lifetimes. When we consider si- 
multaneous heating by both AGN, based on the current out- 
burst, and thermal conduction, the suppression factor nec- 
essary to satisfy the heating requirements within the cooling 
radius is roughly 0.16 which is acceptable. 

Virgo: Even allowing for mass deposition, thermal con- 
duction alone is insufficient for providing the heating re- 
quirements at all radii. The required AGN duty cycle, in 
the absence of thermal conduction, of roughly 10 6 yrs is 
very short compared with predicted AGN lifetimes suggest- 
ing that significantly larger AGN outflows are required at 
certain intervals to prevent a cooling catastrophe. This is 
interesting since Virgo has the smallest heating requirement 
in our sample. Within the cooling radius, the combined ef- 
fect of the current AGN outburst and thermal conduction 
can achieve the required heating with a suppression factor 
of 0.7. 

Perseus: It seems that in the absence of an AGN, ther- 
mal conduction may be insufficient compared to heating re- 
quirements at roughly roughly the central 30 kpc, but may 
be sufficient outside this radius. Without thermal conduc- 
tion, the required AGN duty cycle of roug hly 3 x 10 7 yrs 
is also compatible with predicted AGN lifetimes. The sup- 
pression factor for the combined effect of AGN and thermal 
conduction to provide the required heating is 0.05. 

Hydra: Thermal conduction alone cannot provide the 
necessary heating for maintaining a steady-state. The re- 
quired AGN duty cycle in the absence of thermal conduc- 
tion is long compared with predicted AGN lifetimes which 
reflects the powerful nature of the current outburst. This is 
the most powerful central AGN outflow in this small sam- 
ple. Since the current AGN outburst in Hydra is so power- 
ful, it more than satisfies the heating requirements within 
the cooling radius. Hence we derive a negative value for the 
suppression factor. 

A2597: It appears that thermal conduction on its own 
cannot provide sufficient heating except between roughly 
4-10 kpc. Without thermal conduction, the required AGN 
duty cycle is also compatible with predicted AGN lifetimes. 
A2597 has the largest heating requirements in this sample, 
followed by A478. Both require more than two orders of mag- 
nitude more heating power within the cooling radius than 
Virgo. It is interesting to note that even a combination of the 
current AGN heating rate and full Spitzer thermal conduc- 
tion provide roughly only one fiftieth of the required heating 
rate within the cooling radius. It therefore seems as if the 
only means by which a cooling catastrophe can be averted 
in this system is by an extremely powerful AGN outburst, 
or a cluster-scale merger. 

A1795: On its own, thermal conduction is sufficient to 
match the heating requirements throughout most of the clus- 
ter, although not near 5 kpc. The required AGN duty cycle, 
in the absence of thermal conduction, of roughly 4.3 x 10 7 
yrs is also compatible with predicted AGN lifetimes. For a 
combination of heating by AGN and thermal conduction a 
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Name L^MlO^ergs- 1 ) r/(10 8 yr) ^^/(lO^ergs" 1 ) Lh^AlO^ergs- 1 ) / 



Virgo 


u.uoz_ Q32 


0.078 


13.7 


10.7 


0.72 


Perseus 


o.a_ 1 2 


0.28 


4360 


236 


0.05 


Hydra 


07 i +1-5 
z '- i -2.9 


1.9 


112.8 


231 


-1.8 


A2597 


r: fi +5.2 
°-°-1.6 


0.21 


71.0 


427 


4.8 


A2199 




0.14 


797 


149 


0.16 


A1795 


i2.5j;i 6 3 8 


0.43 


1280 


465 


0.2 


A478 


1 Q+1-0 
J - a -0.3 


0.026 


2720 


1169 


0.1 



Table 5. Assumed average mechanical luminosity per bubble in each cluster based on the observations 
summarized in Birzan et al. (2004) and assuming a duty cycle of 10 s yrs in (column 2), duty cycle of each 
AGN in units of 10 s yrs evaluated at the cooling radius (column 3), rate at which energy is supplied to the 
region within the cooling radius by thermal conduction (column 4), the total required heating rate within the 
cooling radius (column 5) and the suppression factor for combined heating by AGNs and thermal conduction 
(column 6). All quantities are evaluated at the cooling radius. 



suppression factor of roughly 0.2 is required to achieve the 
necessary heating. 

A478: Thermal conduction alone should be sufficient for 
matching the heating requirements. In addition, the required 
AGN duty cycle in the absence of thermal conduction is very 
short compared with predicted AGN lifetimes. The central 
AGN outburst in this system is the least powerful in our 
sample. 

These results suggest that heating from AGN, with duty 
cycles of 10 s yrs, alone is unable to maintain a steady-state, 
but may result in a quasi steady-state if more powerful out- 
bursts are interspersed with the more common lower power 
outbursts. Alternatively, it may be that a combination of 
AGN heating and thermal conduction could provide the nec- 
essary heating, at least within the cooling radius, in all but 
one cluster of our sample: A2597. However, although thermal 
conduction is a possible heating mechanism of the central 
regions of galaxy clusters it is still not clear whether ther- 
mal conduction actually is an efficient heating mechanism 
in galaxy clusters. For examp le, from num erical simulations 
of the Virgo cluster bv lPope et al.l i2005tl . taking in to ac- 
count thermal conduction with different suppression factors, 
the observed temperature profiles are most consistent with 
a suppression factor of between 0-0.1. This is at odds with 
the value of the suppression factor we have derived for most 
clusters in this sample. In all cases the effectiveness of ther- 
mal conduction must be fine-tuned as a function of radius 
in order to allow a steady-state of the cluster gas. 

It is worth noting that galaxy clusters also need to be 
heated at radii outside the cooling radius if a time-averaged 
steady-state is to be maintained. Two possible mechanisms 
for this are sound waves excited by buoyantly rises bubbles 
and comparatively rare, but very powerful, AGN outbursts. 

It may also be possible to predict which clusters are 
more likely to experience more powerful AGN outbursts in 
the near future. For example, if a cooling flow has formed, 
the fractional decrement in gas temperature between the 
cooling radius and the cluster centre is likely to be larger 
than in a cluster which has recently been heated by a large 



AGN outburst. Then, if the duty cycle that we have calcu- 
lated based on the current heating rate, required to maintain 
a steady state, is much less than 10 8 yrs, then we may predict 
that a powerful outburst is iminent. 

Based on the above criteria it seems that in our sample 
A478 and A2597 are the clusters most likely to experience 
a powerful AGN outburst in the near future. Virgo appears 
also to be a candidate for an iminent ouburst, however, the 
relatively small temperature difference between the cluster 
center and the cooling radius indicates that it may recently 
have experienced sufficient heating. 

Our results for the required heating rates indicate that, 
for our initial estimate of the AGN duty cycles, if buoyant 
plasma bubbles are the only heating mechanism by which 
clusters are heated, then, of our small sample, only the AGN 
at the centre of Hydra is currently injecting sufficient energy 
to heat the gas within several multiples of the cooling radius. 
For the remaining six clusters, the energy available in the 
form of plasma bubbles is only sufficient to provide the re- 
quired heating within a relatively small fraction (roughly 
0.1-0.5) of the cooling radius. This significant energy deficit 
implies that these clusters may not be in a steady-state, or 
that the AGNS are active far more frequently at the centres 
than outside of clusters. 

A possible scenario which is consistent with our results 
is that AGNs in general produce only relatively small or 
moderately powerful outbursts, such as we observe in the 
majority of clusters, which delay the occurence of a cool- 
ing catastrophe by a small amount each time. These small 
heating events may occur too infrequently to prevent the 
formation of a cooling flow and the temperature gradient 
grows with time. Eventually sufficient material is accreted 
by the central galaxy to initiate an extremely energetic out- 
burst such as those seen in Cygnus-A and Hydra- A. This 
would also allow the heating of material at radii which are 
hard to reach with lower power outbursts. In this way the 
moderately powerful AGN outflows observed in most clus- 
ters only delay more energetic events. Such a scenario for 
clusters, small-scale AGN heating interpersed occasionally 
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with more energetic events, evolution is complex and would 
have to be modelled numerically. 
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